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ABSTRACT 

We present a study of stellar population of Lyman Alpha Emitters (LAEs) at z = 4.86 in the Great Obser- 
vatories Origins Deep Survey (GOODS) North field and its flanking field. The LAEs are selected based on 
optical narrowband (NB71 1) and broadband (V, I c , and z') observations by Suprime-Cam attached on Subaru 
Telescope. With the publicly available IRAC data in GOODS-N and further IRAC observations in the flanking 
fields, we select five LAEs which are not contaminated by neighboring objects in IRAC images and construct 
their observed spectral energy distributions (SEDs) with I c , z', IRAC 3.6/im, and 4.5^m band photometry. The 
SEDs cover the rest-frame UV to optical wavelengths. We derive stellar masses, ages, color excesses, and star 
formation rates of five LAEs using SED fitting method. Assuming the constant star formation history, we find 
that the stellar masses range from 10 8 to 10 10 M Q with the median value of 2.5 x 10 9 M Q . The derived ages 
range from very young ages (7.4 Myr) to 437 Myr with a median age of 25 Myr. The color excess E(B-V) 
are between 0.1-0.4 mag. Star formation rates are 55-209 M©yr _1 . A comparison of the stellar populations 
is made between three LAEs and 88 LBGs selected at the same redshift, in the same observed field, and down 
to the same limit of the rest-frame UV luminosity. These three LAEs are the brightest and reddest samples 
among the whole LAE samples at z = 4.86. The LAEs distribute at the relatively faint part of UV-luminosity 
distribution of LBGs. Deriving the stellar properties of the LBGs by fitting their SEDs with the same model 
ensures that model difference does not affect the comparison. It is found that the stellar properties of the LAEs 
lie on distributions of those of LBGs. On average, the LAEs show less dust extinction, and lower star formation 
rates than LBGs, while the stellar mass of LAEs nearly lies in the middle part of the mass distribution of LBGs. 
However, the stellar properties of LAEs and LBGs are similar at the fixed UV or optical luminosity. We also 
examine the relations between the output properties from the SED fitting and the rest-frame Lya equivalent 
width, but cannot find any significant correlation. 
Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift 



1. INTRODUCTION 

There are two popular techniques for isolating galaxies 
at high redshift (z > 3): broadband selection (i.e., Lyman 
Break technique) and narrowband selection. Lyman Break 
Galaxies or LBGs are selected by the Lyman Break tech- 
nique by taking the advantage of the spectral discontinu- 
ity due to the neutral hydrogen attenuation in the inter- 
galactic medium at the rest-frame wavelength shorter than 
912/1216 A. Via this method, a large number of high red 
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shift galaxies are studied (e.g., Steidel et al. 1996 
2003; Giavalisco et al. 1996, 1998; Lowenthal et al 
Pettin ietal 111998b IS haplev et al. 20 011 120031: llwata et al I 
2003b 12007b Eap ovich et al. 2004; Bou wens et al. I 2004. 
2008; Reddv et al. 2005, 2006, 2008; Sawicki & Thompson 
2005J 12006b " lYos hida et al71 l2006l) . The other method is to 
select galaxies with strong Lya emission lines which fall 
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into a narrowband filter. This method is useful for select- 
ing high-redshift galaxies that have a strong Lya emission 
line. These ob jects are called Lyman A lpha Emitters (here- 
after LAEs). Partridge & Peebles ( 1967) proposed that pri- 
mordial galaxies in the early stage of their formation should 
show a strong Lya emission line. LAEs are thus expected to 
be young galaxies with low metallicity. Many surveys have 
been made to seach for galaxies with a strong Lya emis- 
sion at various redshifts ranging from 2.1 to 6.6 or even more 
(e.g. [Hu&McMahon 1996; Cowie & Hu 1998 : iHu et al. I 
1998, 2004; Rhoa ds et all 12000b 12004b Kudritz ki et al. 
2000; Rhoads & Malhotra 2001; Malhotra & Rhoads 2002 
Fvnb o et al. 1 12001b lOuchi et al 1 120031 12008b iFuiita et al. 
2003b Shimasa kuet al. 1 120031 120041 12006b iKodaira et al 
2003b lAiiki et al 1 120031 [2004, 2006: Taniguchi et al. 2005; 
Venemans ~al. I 12002b l200l live et al. I 120 06: Otaetal. 
2007; Nil sson et al. ll2009blGuaita et al. II2009IX 

In order to examine the evolutionary stage of LAEs 
and to reveal what kind of galaxy an LAE is, various 
properties of LAEs have been studied such as luminosity 
funct io ns (e.g.,[O uchi et al. 2003, 2008; Maihotr a"& Rhoads I 
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2004b ISh imasaku etal. I 12006b ,„,,„„„,,„,,>.,„. 
Gronwall et al. 112007b lOta et al. Il2007l) or clustering 



ties ( e. g.. lOuchi et al" 1 12003b ' Shimasaku et al. 2003, 2004, 
120061: Kashikawa et al. |20Q6i IMurayama et al. 2007). 
Revealing the stellar population of LAEs is one of crucial 
studies for understanding their physical nature. In order to 
do that, spectral energy distribution (SED) of a galaxy is 
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compared with stellar population synthesis models produced 
by varying the ages, metallicities, amounts of dust extinc- 
tions, SFRs, etc. The stellar population of a galaxy can be 
constrained from the best fit model. This method is known 
as SED fitting method. Recent studies show that LAEs are 
in wide ranges of ages (lMyr - lGyr) and stellar masses 
(10 6 - 10 10 M ? ) (e.g.. IGawiser et al 1 120 06. 200 71: lLai et al 

2007, 2008; Nilssonetal. 2007; Finkelstein et al. 2007 

2008, 2009; Pirzkal et al. I |2007|) . Among these studies, 
stacking analysis shows that LAEs at z ~ 3 are free from 
dust or show modes t dust extin ction (E(B-V) < 0.03 mag) 
(e.g. IGawiser et al. 1 2006. 2061; lLai et al~ll200l lOno et al. I 
2009), whereas SED fitting of an individual LAE at higher 
redshifts sugge sts that some of LAEs show significan t dust 
extinction (e.g jLai et al. Il2007t iFlnkelstein et al. 1120091) . 

It is important to investigate the connection between LAEs 
and other galaxy populations selected by different method, 
i.e., LBGs that are selected based on UV continuum. Since 
these two methods suffer from different biases, the selected 
galaxies partial ly overlap and the relationship between them 
is not clear. Shapley et al. I d2001l) divided LBGs at z ~ 3 into 
2 subsamples according to their ages obtained by SED fitting: 
"old subsample" with ages about 1 Gyr and "young subsam- 
ple" with ages less than 35 Myr. They found that the old sub- 
sample shows a strong Lya emission line, while young sub- 
sample do not. In other words, they found that LBGs with 
Lya emission line are older than those without Lya emission 
line. Confirming this statement, the recent wor k with a larger 
sample at the same redshift by iKornei et al. I ((2009) shows 
that objects with rest-frame Lya equivalent width larger than 
20 A seem to be older, lower in star formati on rate, and less 
dusty than those without Lya emission line. Pent ericci et al. I 
(2007), in contrast, found that at z ~ 4 LBGs with Lya emis- 
sion are less massive and younger than those with no line. 
By considering the results from the SED analysis mentioned 
above, there may be an evolutio n of difference between LAE s 
and LBGs with redshifts (e.g., Shimiz u~& Umemurall2010l) . 
Although there are some SED studies at higher redshift (z > 5) 
(iPirzkal let al. I l2007t lLai et al. 1 120071) . the relationship be- 
tween LAEs and LBGs is still unknown. In this paper, we 
study the stellar populations of LAEs at z = 4.86 from their 
rest-frame UV-to-optical spectral energy distributions. Down 
to the same UV luminosity, the derived stellar properties o f 
LAEs are compared to those of LBGs bv lYabe et al. I (2009). 
They selected LBGs at the same redshift range in the same 
field and they used the same fitting code and synthesis spectral 
model assumptions. This direct comparison between LAEs 
and LBGs is expected to reveal their difference and to know 
what kind of properties make a galaxy an LAE. 

Data sources and photometry are described in section 
12 Section [3] explains the selection criteria, and the LAE 
candidates. The observed SEDs are constructed for each 
LAE candidate and compared to the models in section[5] The 
fitting results are shown in Section [6] Section [7] is for com- 
parisons between LAEs and LBGs. The summary is given 
in section [8] Through out this paper, we use AB magnitude 
system (Oke & Gunn 1983) and adopt a cosmology with 
parameters of il m = 0.3, Oa = 0.7, and Hq = 70 km s" 1 Mpc -1 . 

2. DATA SOURCES AND PHOTOMETRY 

2.1. Optical data 
2.1.1. Obser\'ations and data reduction 
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Figure 1. Observed fields. The solid line indicates the effective field ob- 
served with Suprime-Cam, whereas the long and short dashed lines show the 
GOODS-N and GOODS-FF fields observed with IRAC. Hubble Deep Field 
North (HDF-N) is shown as a reference at the center of the figure. 

Optical data we re obtained with Suprime-cam 
(Mivazaki et al. I [20021) attached on the Subaru telescope 
dive et al. 1 12004]) . The observed fi eld is toward th e Hub- 
ble Deep Field-North (HDF-N; iWilliams et all [19961) 
[RA(2000)=12 h 36 m 49. s 4, Dec(2000)=+62°12'58"] and is 
illustrated by solid line in Figure Q] Pixel scale of the CCD 
was 0."20. We used NBlll [X eff = 7126, FWHM = 73] 
and three broadband filters: V, I c , and z' to select z = 4.86 
LAEs. The effective redshift interval calculated from FWHM 
of NBlll is 4.83 < z < 4.89. Transmission curves, which 
include the responses of CCD, prime focus corrector (PFC), 
mirror, and airmass, of all filters used in this work are shown 
in Figure The observation with the NB7 1 1 filter was 
made on March 16, 2005. Images were taken with dithering 
of ~ 80" and the exposure time of 1200 seconds for each 
frame. With 23 exposures, we covered ~ 750 arcmin field 
of view. The total integration time was 7.7 hours. The 
weather condition was not so excellent. Seeing size during 
the observation was averagely l."5. 

Data reduction was made by using SDFRED version 1.4.1, 
the software developed for Supr ime-Cam data reduction 
dYagi et al. Il2002t lOuchi et al. 1120041) . After the bias subtrac- 
tion, flat fielding, and distortion correction, image mosaicking 
was made by identifying 30- 170 non-saturated stars common 
in several object frames. The frame alignment and correction 
for flux/count and full width at half maximum (FWHM) were 
made based on these star data. A resulting FWHM of point 
sources in the mosaicked image is ~ l."7. 

Astrometry was made based on USNO-A2 catalog by fit- 
ting with the third-order polynomial coefficients to ~ 1 , 000 
stars identified in the mosaicked image. The positional error 
is about 0.2" rms over the image. The magnitude zero point 
was derived based on the imaging data of two spectrophoto- 
metric standard stars (G191-B2B and HZ44) taken during the 
same observing night. As an independent check for photome- 
try in NBlll images, we derived NBlll magnitudes of stars 
by interpolating their V, I c , and z' magnitudes. The derived 
magnitude zero points are in agreement within ~ 0.18 mag. 
A 3cr limiting magnitude at 2. "5 diameter aperture is 26.16 
mag for the NB1 1 1 image. 

The observations and data reductions of broadband images 
are described in detail by llwata et aLI (120071) . The images 
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Figure 2. Transmission curves of V, NB711, I c ,andz' bandpasses. CCD re- 
sponse, transmission of PFC, mirror reflectivity, and airmass (sec z = 1 .2) are 
included. Also shown is a model spectrum in the unit of fijy (solid line) of a 
star forming galaxy at z = 4.86 including Lya emission line with a rest-frame 
equivalent width of 10A. The model is constructed by Bruzual & Chariot 
(2003) stellar synthesis code. A constant star formation history and Salpeter 
IMF are assumed at an age of 12 Myr and at star formation rate of 1 MQyr~' . 
The attenuation by intergalactic medium is applied with the prescription by 
Madau(1995). 

were taken with V, I c , and z' filters in February 2001 . The typ- 
ical seeing during the observations was ~ l."l. In addition, 
llwata et aLl (|2007l) also used the imaging data taken from 
February t o April 20 01 by the University of Hawaii (UH) 
group (Ca pak et al. 12 004). who used the same instrument and 
filters. The final effective survey area is 508.5 arcmin 2 . The 
limiting magnitudes of V, I c , and z'-band images are 28.1, 
26.8, and 26.6 mag, respectively (3a at l."6 diameter aper- 
ture). The broadband images were all degraded to a PSF size 
of l."7 to match that of the NBlll image. The 3a limit- 
ing magnitudes of final PSF matched images at 2.5" diameter 
aperture are 27.3, 26. 1, and 25.9 mag for V, I c , and z 1 images, 
respectively. 

2.1.2. Photometric catalog 

The photometric catalogs of the optical images were made 
by using SExtractor version 2.5.2 (Bertin & Arnouts 1996). 
The positions of objects were extracted from the NBlll im- 
age and the photometry was then made at the extracted posi- 
tions in all optical images via dual-image mode. In order to 
use the dual-image mode, we first have to make a position 
registration between the broadband images and the NBlll 
image. Position mapping between the images was made by 
using geomap in IRAH3 based on positions of the bright but 
non-saturated stars detected in the images. Then the registra- 
tion of the images was performed by gregister. The accuracy 
in making a position registration is within 0.2 pixels or 0."04 
rms. After image registration, we subtracted GLOBAL back- 
ground from each registered image individualljQ and homog- 

7 IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As- 
tronomy, Inc., under cooperative agreement with the National Science Foun- 
dation. 

8 Based on our simulation with artificial objects, the processes of subtract- 
ing background before smoothing the images and re-subtracting the LOCAL 
background in SExtractor give the better photometry. 



enized the broadband images so that they all have the same 
seeing size as the NBlll image within 0."01 accuracy by 
gauss. The dual-image mode was then performed by using 
the following parameters. The images were filtered with a 
default convolution kernel (default. conv). The LOCAL back- 
ground was estimated using 64 x 64 pixel background mesh 
with 3x3 median filtering. An object was detected with a 
minimum of 5 connected pixels above 2. Oct minimum thresh- 
old. These parameters are found to maximize the number of 
detected objects in the NBlll image and minimize the de- 
tected numbers of objects in the negative NB1 1 1 images. 

In SED fitting process, the total photometry of the objects 
is necessary. The total magnitudes were obtained by applying 
aperture corrections to 2. "5 diameter aperture magnitudes. 
We examined the best value of the aperture size and chose 
2. "5, which maximizes the signal-to-noise ratio (S/N). The 
aperture corrections were determined from the Monte Carlo 
simulations where artificial objects with PSF shape corre- 
sponding to the seeing size of l."7 were put into the original 
image and were then detected by using the same SExtractor 
parameters. In the simulations, our objects were implicitly as- 
sumed to be point sources since their apparent size in the im- 
age is comparable to that of PSF. The PSF was made by stack- 
ing images of the point spread function that shows the stel- 
larity index larger than 0.98 and has no nearby objects. The 
aperture correction factors are -0.30 ± 0.01 and -0.33 ± 0.01 
mag for I c and z' images, respectively, where the errors are 
the PSF uncertainties. Errors of the total magnitudes are the 
combination of la error in making photometry of simulated 
objects and the uncertainties of correction factors. Because 
we do not use V and NBlll images in SED fitting process, 
we do not compute the aperture corrections for them. 

2.2. Mid-infrared data 

Mid-infrared images are obtained from deep observations 
with Infrared Array Camera (IRAC) on Spitzer Space Tele- 
scope (SST). We used the publicly available mid-infrared 
data in the Great Observator ies Origins Deep Survey North 
field jDickinsonetal.1120031 hereafter GOODS-N) centered 
at 12 h 36 m 54. s 87, +62°14'19."2 (J2000) provided by Spitzer 
Space Telescope Legacy Science prograrrfj The field covers 
an area of approximately 10' x 16' or ~ 160 arcmin 2 as shown 
in Figure Q] In this work, we used the First Data Release 
(DR1) and Second Data Release (DR2) of IRAC data consist- 
ing of imaging data in 3.6, 4.5, 5.8, and 8.0 /im bandpasses. 
Pixel scale of all images after drizzled is 0."60 pixel -1 . Mean 
FWHM of PSF in IRAC GOODS-N images is -l."7. 3a lim- 
iting magnitudes of IRAC 3.6, 4.5, 5.8, and 8.0 /im images at 
2. "4 diameter aperture are 25.9, 25.8, 23.7, and 23.6 mag, 
respectively. 

In addition to the very deep GOODS-N data, data in the 
flanking field (hereafter GOODS-FF) of GOODS-N are ob- 
tained to cover a part of the Subaru observation area. The 
GOODS-FF observations were carried out with IRAC cam- 
era in December 2005 and June 2006 (the Spitzer GO pro- 
gram (GO-20218); PI = Ikuru Iwata). The observations were 
made by pointing at 5 different positions on the sky around the 
GOODS-N field. The total exposure time for each passband is 
4,000 seconds. Basic Calibrated Data (BCD) processed by the 
pipeline of Spitzer Science Center (SSC) were used. The im- 
ages were drizzled using MOPEX software provided by SSC 

' http: //ssc. spitzer.caltech.edu/legacy/goodshistory.html 
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and combined together to get an improved mosaic image cov- 
ering the whole observed area of ~ 300 arcmin 2 . The mosaic 
image has a final pixel scale of 0."61 pixel -1 and a PSF size 
(FWHM) of l."7. The 3a magnitude limits at 2."4 diame- 
ter aperture are 25.0, 24.6, 22.1, and 22.3 for IRAC 3.6, 4.5, 
5.8, and 8.0/im bands, respectively. Combining the infrared 
data obtained in GOODS-N and GOODS-FF fields, we have 
the effective area of ^400 arcmin 2 covering ^80% of Subaru 
images as seen in Figure Q] 

The photometry of IRAC was made based on the positions 
of LAEs at z = 4.86, which are selected from the criteria de- 
scribed in the next section. Since the dimensions and the pixel 
scale of the narrow-band image and IRAC images are differ- 
ent, we could not use dual-image mode in SExtractor to de- 
termine aperture photometry in IRAC images. Phot provided 
by IRAF was used instead to make aperture photometry at 
2. "4 diameter, which are examined to be the best aperture size 
since it maximizes the S/N ratio. Because the positions in the 
NB1 1 1 image are used as the reference, we made a registra- 
tion of IRAC images before making IRAC photometry of the 
objects. The position alignment was made globally for each 
IRAC image by using positions of point sources detected in 
both NBlll and IRAC images. The estimated errors in de- 
termining the position shifts are 0."1 (0.17 pixels) for IRAC 
images in GOODS-N field and 0."2 (0.33 pixels) for those 
in GOODS-FF. Since the error of 0."2 may not be negligi- 
ble, we re-made the photometry of objects at the positions 
separated from the corrected positions by 0."2 and found the 
magnitude difference to be less than 0.2 mag. This magnitude 
difference gives less than 3% difference in the stellar masses 
derived from SED fitting. 

Aperture corrections for IRAC images were determined by 
the following process. The artificial objects with IRAC PSF 
are put into the IRAC images. The PSF was made from 
images of objects identified to be a stellar object in Sub- 
aru images. For the GOODS-N field, the correction factors 
are -0.62 mag and -0.72 mag in 3.6/mi and 4.5/im bands 
respectively. The correction uncertainties due to the uncer- 
tainties in PSFs are 5% and 6% in 3.6/mi and 4.5/im bands 
respectively. For GOODS-FF images, they are respectively 
-0.70 mag and -0.73 mag with 6% and 7% uncertainties in 
3.6/im and 4.5/im bands. Errors in IRAC photometry are la 
standard deviations of sky background and the uncertainties 
of correction factors. Since 5.8/im and 8.0/im images show 
too low signal-to-noise (S/N) ratio to place a useful upper 
limit in SED fitting process, we do not use them in further 
analysis. 

3. SAMPLE SELECTION 

29,675 objects were detected with NBlll < 26.0 mag 
(3.5cr limiting magnitude). Figure [3] shows I c — NBlll color 
versus NB71 1 magnitudes of all objects (black dots). We se- 
lected NBll 1-excess objects from the following selection cri- 
teria: 



I c -NBll\> 0.64 and 

I c -NBlll>3a Ic . NBln , 



(1) 

(2) 



where 3ai c -MBin is a 3a error in measuring I c -NBlll color 
for a source with a color of I c -NBll I = -0.28 mag, the aver- 
age color of all objects with NB7 1 1 magnitudes brighter than 
26.0 mag. Note that the magnitudes and colors of objects 
used in the sample selection are the values at 2. "5 diameter 
aperture, except for those noted otherwise. The first criterion 



2.5 

2 

1.5 

1 

0.5 



-0.5 

-1 



1 ' ' ' 1 ' ' ' ' 1 ' ' \ 
z' >= 26.5 * 


,1111,111 


' ; 


z' < 26.5 • 




: 




• \ 


'■ 


: 


' • *\ 


/ ■ 


- : 


\ i 


7 , 


\ - 

-— -- -A: 


- 


: / v ., ...... ; 


. i . . . . i ;■ . . . 1-. .' . , - 



21 



22 



23 



24 
NB711 



25 



26 



27 



Figure 3. I c -NBlll color versus NBlll 2. "5 aperture magnitude. LAEs 
with z t ' otal < 26.5 mag are shown with rilled circles; rilled triangles represent 
those with z' . > 26.5 mag. Arrows show the lower limits of I c -NBlll 
taken at 2<r I c limiting magnitude. Note that we used 2a I c limiting magni- 
tude as a lower limit for objects whose l c magnitudes are fainter than that. 
A diagonal black line indicates 2a limiting magnitude of the I c image. Red 
curves indicate the distributions of 3a errors in measuring I c — NBlll color 
for a source with a color of I c — NBlll = —0.28 mag, the average color of 
all objects brighter than NBlll = 26.0 mag. The horizontal and vertical ma- 
genta lines show the NBll 1-excess criterion and NBlll = 26.0 mag limit, 
respectively. 

corresponds to a rest-frame equivalent width (EW rest ) of 10 
A assuming a flat continuum (f v = constant). The EW cut 
of 10 A used in this paper is lower than that commonly used 
in other LAE studies (20 A ). Other studies (e.g., Ouc hi et al. I 
120031: iNilsson et al. 1120091) compute the EW by using the av- 
erage continuum flux density from the continuum both blue- 
ward and redward of the Lya line, while the continuum in 
this work is extrapolated from the continuum redward of the 
Lya line. Because the blueward continuum suffers from the 
attenuation by neutral hydrogen in the intergalactic medium 
(IGM) at high redshift, the continuum estimated by the lat- 
ter way is expected to be higher than that in the usual case. 
Consequently, the EW computed in this paper is smaller than 
that from other studies, even though the object has the same 
redward continuum. The Lya EWs computed in the usual 
way are also shown in Table |2 most of the LAEs compared 
to LBGs show the EWs larger than 20A when we adopt the 
commonly used way. 

By the above criteria, 667 objects were selected. Among 
them, there are low-redshift interlopers with strong emission 
lines such as Ha, UB, [OIII]A5007, or [OII]A3727. To rule 
them out, we first divided objects selected by the above cri- 
teria into two groups according to their brightness in z' band. 
For candidates with z t ' otal < 26.5 mag, we adopted the follow- 
ing criteria, which were used by Iwata et al. (2007; see their 
figure 2) to select LBGs at z ~ 5: 



Ztotal< 26 - 5 , 

V-I c > 1.55 and 
V-/ c >7.0(/ c -z') + 0.15, 



(3) 

(4) 
(5) 



where z t ' oml is the total magnitude in z' band and colors are de- 
rived from 2."5-diameter aperture magnitudes. These criteria 
have been confirmed by spectroscopic follow-up observations 
dAndo et al. 1 12004. 2007; iKaiinoetall 120091) . For objects 
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Figure 4. Montage images of group I LAEs. Size of each image is 10" X 10" 
with an LAE at the center of the image. North is at the top; east is to the 
left. IDs of objects are shown on the left of the figure. Note that the Subaru 
broadband images displayed here (V, I c , and z') are those before smoothing. 
An asterisk indicates the LAE with the spectroscopic redshift. 
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Figure 5. Same as Figuref4]but for group II LAEs. 
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Figure 6. Same as Figuref4]but for group III LAE 



with Zj otal > 26.5 mag, we required non-detection in V-band 
magnitude (V > la limiting magnitude). As a result, 24 ob- 
jects meet all criteria simultaneously down to NB1W = 26.0 
mag as illustrated in Figure [3] Eleven of them show z' mag- 
nitudes brighter than 26.5 ma g. Note t hat zL ,»\ here is slightly 
different from z„ agauto bv llwata et al. I (120071) (up to ±0.3 mag 
at z 1 = 26.5 mag) because the broadband images used in this 
work were smoothed to match the NB711 seeing size before 
making the photometry. 

In order to study the properties of LAEs by SED fitting, 
the rest-frame optical photometry (i.e., IRAC photometry) is 
needed. 16 out of 24 LAEs are in IRAC images (hereafter 
IRAC sample). Since the angular resolution of IRAC im- 
ages is not excellent, some LAEs are contaminated by nearby 
objects. In order to examine whether or not neighboring 
objects affect photometry of LAEs, we simultaneously fit- 
ted all neighborin g objects within 10" radius with GALFIT 
dPeng e t al. 2002) and subtracted them from the original im- 
age. In the GALFIT process, we masked the LAEs out so 
that the accurate sky estimation can be obtained. Then we 
made aperture photometry of the LAEs from the residual im- 
ages. We only selected LAEs of which magnitude difference 
between the original images and the residual images are less 
than 0.1 mag as isolated LAEs. Finally, we have 12 isolated 
objects. We divided these 12 LAEs into four groups based 
on their brightness in rest-frame UV and optical wavelengths. 
Group I is for two LAEs which are detected above 2a limit- 
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Figure 7. Same as Figure|4]but for group IV LAEs. 

ing magnitudes in all bands used in SED fitting process (i.e., 
I c , z', IRAC 3.6^m and 4.5^m). One of them, #35915, is 
spectroscopically confirmed to be at z = 4.82. The emission 
line does not show a doublet profile in R ~ 2000 spectrum ob- 
tained by GMOS on Gemini telescope (Ouchi et al. in prep). 
With the given resolution, this can rule out the possibility that 
the object is an [Oil] emitter. Moreover, there seems to be a 
continuum break at the wavelength shortward of the line. Two 
LAEs are classified as group II which is an LAE detected in 
two bands: one in Subaru optical and one in IRAC bands. 
One LAE is classified as group III which is not detected in 
I c and zf bands but detected in both 3.6/im and 4.5/im bands. 
Group IV is for 7 LAEs that are not detected above 2a limit- 
ing magnitudes in three or more bands. In SED fitting process, 
it is necessary that an LAE should be detected above 2a mag- 
nitude limit in more than two bands; therefore, we can use 
only group I— III LAEs. Hereafter we call group I— III LAEs 
SED fitting sample. The montage images of LAEs are shown 
in Figure [4] [5] and |7j for group I, II, III and IV, respec- 
tivelyH Aperture photometry of the isolated LAEs is summa- 
rized in Table Q] Neighboring objects are seen around some 
of group I-II LAEs, but note that the neighboring objects do 
not affect the photometry more than 0.1 mag. It is noteworthy 
that the object #23098 in group III is undetected in all optical 
bands except for the NB711 image, while a bright object is 
seen close to #23098 with a slight offset in IRAC 3.6/im and 
4.5yum images. Our aperture photometry includes the light 
from this object. However, we cannot rule out the possibility 
that the object is a foreground red object. When we subtract 
the object, the aperture magnitudes in IRAC bands are fainter 
than 2a limiting magnitudes. Thus it should be kept in mind 
that the SED of this object may be contaminated. 

4. SELECTION BIAS 

10 The adjacent object of #22962, which is very bright in NB711 band, is a 
separate object and is likely to be a low-redshift emitter, because it is detected 
in V band and does not satisfy the V -I c > 1.55 criterion. 
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Table 1 

Aperture photometry of isolated LAEs 



ID 


Field 


ya,b 


NB711 a 


/ c a,c 


_/a,c 


3.6Mm a - d 


4.5/^m a - d 






(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


Group I 


30702 


GOODS-N 


28.37 ±0.69 


25.46±0.18 


26.25 ±0.36 


26.12±0.36 


24.97 ±0.14 


25.36 ±0.22 


35915 


GOODS-N 


> 28.51 


25.06±0.12 


25.79 ±0.25 


25.61 ±0.24 


25.06±0.16 


25.60 ±0.27 


Group II 


12766 


GOODS-FF 


> 28.51 


25.81 ±0.23 


> 26.50 


26.33 ±0.42 


25.12±0.35 


> 25.05 


22962 


GOODS-N 


> 28.51 


25.62 ±0.20 


26.26 ±0.37 


> 26.38 


26.06 ±0.35 


> 26.24 


Group III 


23098 


GOODS-FF 


> 28.51 


25.85 ±0.24 


> 26.50 


> 26.38 


24.42 ±0.20 


24.21 ±0.23 


Group IV 


08550 


GOODS-FF 


> 28.51 


24.39 ±0.07 


26. 14 ±0.33 


> 26.38 


> 25.42 


> 25.05 


09296 


GOODS-FF 


> 28.51 


25.83 ±0.24 


> 26.50 


> 26.38 


> 25.42 


> 25.05 


11583 


GOODS-FF 


> 28.51 


25.57±0.19 


> 26.50 


> 26.38 


> 25.42 


> 25.05 


13190 


GOODS-FF 


> 28.51 


25.78 ±0.23 


> 26.50 


> 26.38 


> 25.42 


> 25.05 


14814 


GOODS-FF 


> 28.51 


25.84 ±0.24 


> 26.50 


> 26.38 


> 25.42 


> 25.05 


15089 


GOODS-FF 


> 28.51 


25.27±0.15 


> 26.50 


26.30 ±0.42 


> 25.42 


> 25.05 


27510 


GOODS-N 


> 28.51 


25.56±0.19 


> 26.50 


> 26.38 


> 25.42 


> 25.05 



a Errors are lcr values. 

b Upper limits are \a values at 2. "5 diameter aperture. 
c Upper limits are 2<r values at 2. "5 diameter aperture. 
d Upper limits are la values at 2. "4 diameter aperture. 

In this section, we check the sample bias in two points 
of views: bias on selecting the isolated LAEs and bias on 
choosing the SED fitting sample. Figure |8ja) shows the dis- 
tributions of rest-frame Lya equivalent widths of all LAEs, 
IRAC sample, isolated LAEs (group I-IV), SED fitting sam- 
ple (group I— III), and the compared LAE sample (see section 
IT). The rest-frame Lya EWs are calculated by assuming that 
all LAEs are at z = 4.86 and the Lya emission line falls into 
the center of NBll 1 band (Tabled. Note that the values may 
be underestimated if the LAEs are not exactly at z = 4.86, i.e., 
the Lya emission line does not fall into the center of NB711 
band. In addition, the EWs of 20 objects (of all 24 LAEs) 
are more underestimated due to the 2a upper limit on I c mag- 
nitudes. For the spectroscopically confirmed LAE (#35915), 
the Lya emission line is in the NB1W band, but is out of 
NB7 1 1 FWHM range. By considering the transmission curve 
of the NB1 1 1 filter, the corrected EW is estimated to be 40 
A (Table |2]i and is used hereafter. As a cross-check, we also 
examine the EW from the spectrum and find that it is consis- 
tent with that obtained from the imaging. The EWs of LAEs 
distribute from 10 A (the selection limit) to 55 A. Figures[HJb) 
and (c) show the distributions of z' magnitudes and I c —z' col- 
ors, respectively. Leftmost bins in Figures [8] (b) and (c) rep- 
resent LAEs with z' magnitudes fainter than 26.5 mag and 
with magnitudes fainter than 2a limiting magnitudes in both 
I c and z' bands, respectively. The large number of LAEs in 
the leftmost bins of both figures imply that most of LAEs have 
faint UV continuum. As seen in Figures[8](a)-(c), the isolated 
LAEs (red open histogram) occupy the same range of distri- 
butions as all sample at z = 4.86 (orange histogram). In order 
to check whether or not selecting isolated sources in IRAC 
images can cause any bias on properties of the whole LAEs, 
we applied the Kolmogorov-Smirnov test (KS test) to the dis- 
tributions of properties of all LAEs at z = 4.86 and the isolated 
LAEs. We could not reject the null hypothesis that the rest- 
frame Lya EW, z 1 magnitude, and l c -z' color distributions 
of these two samples are drawn from the same populations at 
more than 95% confidence levels. However, we could not use 
all isolated sample in SED fitting because of their faintness 



either in rest-frame UV or optical wavelengths. To investi- 
gate if SED fitting in the following section can represent all 
populations of LAEs at this redshift, we have to compare the 
SED fitting sample (group I— III, green histogram) to all LAEs 
(orange histogram). It is seen in Figur4H](a) that group I— III 
LAEs represent both high- and low-EW LAEs. However, we 
can see from Figures [8] (b) and (c) that the SED fitting sam- 
ple are biased toward the bright UV luminosity and red I c —z' 
colors as we neglect the leftmost bins where magnitudes are 
unreliable. In addition, Table[T]shows that SED fitting sample 
is relatively bright in IRAC bands as compared to all isolated 
LAEs implying the brighter luminosity in rest-frame optical 
wavelength. Thus we conclude that the SED fitting and the re- 
sults hereafter are for the LAEs which have the brighter rest- 
frame UV and optical magnitudes, and the relatively redder 
UV colors. 

5 . STELLAR POPULATION SYNTHESIS MODEL AND SED FITTING 

In this paper, we intend to compare stellar populations of 
the LAEs to those of LBGs at the same redshift by Yabe et al. 
(2009). In order to make a fair comparison, we us ed the same 
stellar population synthesis model as that used in Yabe et al. 
(2009). The model SEDs were obtained by Bruzual & Chariot 
synthesis code (2003, hereafter BC03). We used Padova 1994 
evolutionary track as recommended by BC03. Salpeter(1955) 
initial mass function (IMF) with lower and upper mass cutoffs 
of 0.1 and lOOM is assumed. We made models by fixing the 
metallicity at 0. 2Zd n \ and assuming the constant star forma- 
tion history. BC03 uses quasi-logarithmic 221 time steps from 
0. 1 Myr to 20 Gyr. Time steps were adopted to 5 1 logarithmic 
steps both to reduce the calculation time and to avoid dealing 
with an unequally spaced scale of the original 221 models. 
The age of the Universe at z ~ 5 is ~ 1 .2 Gyr. However, as a 
cross-check on the fits, we allowed the age up to the oldest one 
available in BC03. The effect of dust attenuation is taken into 

1 ' Fixing the metallicity at the lower abundance (0.005Zq ) does not sig- 
nificantly change the fitting results. The stellar masses derived from the lower 
metallicity model differ from those derived from 0.2Zq model by ~ 10% at 
most. The average differences of age, color excess, and star formation rate 
are ~ 50%, ±0.2 mag, and ~ 5%, respectively. 
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Figure 8. Distributions of rest-frame Lyo EW, z'-band magnitude, and I c —z! color of LAE samples. We divide each figure into 2 panels for display purpose. In 
the top panel, an orange histogram shows the distribution of all 24 LAEs at z = 4.86, while the blue and red histograms are for the IRAC samples and the isolated 
LAEs (group I— IV), respectively. In the bottom panel, red histograms refer to the isolated LAEs (the same as in the top panel). Group I — III LAEs are shown 
in green histograms. Distributions of group I LAEs and #12766 from group II LAEs (the compared LAE sample; see section[7} are shown in grey histograms. 
Arrows shown in panel (a) represent the lower limits in Lya EWs of LAEs; the number of arrows directly corresponds to the number of LAEs with Lya EW 
lower limits. Note that the leftmost bins in panels (b) and (c) indicate, respectively, the LAEs fainter than z! = 26.5 mag and LAEs that are not detected above 2rr 
limiting magnitudes in neither I c nor z! bands. 



account by using Calzetti extinction law (Calzetti et al. 2000) 
by changing the color excess E(B-V) from 0.0 mag to 0.8 
mag with a 0.01 step. Attenuation by the IGM is calculated 
with the prescription by Madau (1995). The model spectra 
were then convolved with the appropriate filter transmission 
curves to give model fluxes. This is exactly the same as done 
in Sawicki & Yee (1998). Except for one spectroscopically 
confirmed LAE (#35915; z = 4.82), the redshift was fixed at 
z = 4.86 under assumption that Lya emission of our LAE 
sample is detected at the center of the NB71 1 bandpass. 

Some of the observed SEDs including one spectro- 
scopically confirmed LAE show a significant excess in 
the 3.6/im band as compared with the magnitude in the 
4.5/imband. This is likely to be due to the Ha emission line, 
which falls into IRAC chl (3.6/im) at this redsh ift. The ex- 
cess i s also seen in some fraction of z ~ 5 LBGs ( Yab e et aL~l 
120091) . The SEP model including Ha emission shows a bet- 
ter fit dCharvet al. 112 005: Finkelstein et al. 2008; Yabe et al. I 
120091) . As adopted by lYabe et al. I (120091) . the spectrum of 
Ha emission is included in the synthesis model spectrum by 
the following process. The luminosity of Ha emission is cal- 
culated from the star formation rate of the model by using the 
relation by Kennicutt (1998)0 The dust extinction to the line 
(Calzetti et al., 2000) is assumed to be the same as that to 
stellar componentQ The Ha flux density is finally put into 
each model SED. The existence of Ha emission line in the 
models makes the fit better without adding a free parameter 
(lYabe et al. 1120091) . 

We use the SEDfit software (Sawicki, in prep.), which 
is an evolved v ersion of the SED-fitting software used in 
Sawicki & Yee ( 1998) and subsequent papers, including the 
z ~ 5 LBG study by Ya be et al. I (120091) . After we produce 
a suite of model fluxes by means described above, the best- 
fitting model is found by means of a maximum likelihood test 
for each object. For objects detected in all bands (i.e., our 
group I objects), this test is the standard \ 2 minimization. For 
objects that are undetected in some of the bands (i.e., group 

12 Adopting this relation may be disputable. Dependences on metallicities 
and differences of extinction to stellar continuum and nebular emission are 
discussed by Yabe et al. (2009; their Appendix B). 



II and III objects), SEDfit employs a modification of the \ 2 
formalism: as is standard in the \ 2 approach, for the filters in 
which the object is detected the maximum likelihood calcu- 
lation considers the likelihood that the detected data deviate 
from the model given the uncertainties; for the undetected fil- 
ters, the calculation adds in the likelihood that a given model 
would/would not have been detected given the upper limit of 
the non-detection. For more details on this technique see Saw- 
icki (in prep.). 

In this paper, we constructed an observed SED of each in- 
dividual object from the photometry in I c , z', IRAC 3.6/im, 
and IRAC 4.5/im bands. As explained in section [2] we do 
not use V-band photometry in the SED fitting process in or- 
der to avoid the uncertainty due to IGM absorption. IRAC 
5.8/im and 8.0/tm photometry is not used because of the low 
S/N ratio of the images. Including the upper limits of the pho- 
tometry in IRAC 5.8/im and 8.0/im bands does not usefully 
improve the fitting. Free parameters in our fitting process are 
age, color excess, and scaling normalization. The SFR is ob- 
tained from the scaling normaliz ation. We then find the stellar 
mass by using the age and SFR (Sawic ki & Yee 1 11998). 

6. RESULTS 

The fitting results are summarized in Table [2] Errors in 
the table are at 68% confidence level and are determined as 
follows. The Monte Carlo realizations for each object are 
performed; we vary the input fluxes within their photomet- 
ric uncertainties, rederive the best-fitting model and repeat it 
100 times. The error on an individual parameter is then 68% 
range of the realizations in that parameter. The best-fitting 
model spectra are shown with the observed SEDs in Figure|9| 

The derived stellar masses are ranging from 10 8 — 10 11 
Mq with the median value of 2.5 x 10 9 M©. A typical error 
on the stellar mass is ~ 0.9 dex. By considering only group I 
LAEs which are detected in all bands used in the SED fitting, 
the stellar mass is in the order of 1O 9 M with the smaller 
typical error of ~ 0.6 dex. A relationship between the stel- 
lar masses and the rest-frame optical absolute magnitudes is 
shown in Figure [TOT a); the optical absolute magnitudes were 
calculated from the [4.5/im] magnitudes by assuming f\ oc A' 3 
with f3 derived from the z 1 - [4.5/im] color of an individual ob- 
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Table 2 

The best fit results 



ID Field 


rest Lya EW a 

(A) 


log[Mass] 
(M Q ) 


log [Age] 

(yr) 


E(B-V) 
(mag) 


log[SFR] 

(M yr-') 


A././ 


<? 


Group I 


30702 GOODS-N 
35915 GOODS-N 


14(51) 
40 (202) 


a 6 i+0 .99 
q ^q+0.72 


7 40+2.M) 
7 5Q+1-87 


0.28$ 

o.2o;«?» 


2 15+ 033 
-1.39 

2 03+ 032 


0.02 
0.08 


1.32 
1.16 


Group II 


12766 GOODS-FF 
22962 GOODS-N 


> 11(> 44) 

10(39) 


9 29+ 1 - 59 
8.60*f 


fi 07+2J0 
6.S/_ ag3 


0.30 + " 3 4 
14+kio 


2 - 32+ l'66 
1.74+° : ff 


44.90 
46.63 


< 1.44 

2.27 


Group III 


23098 GOODS-FF 


> 10(> 11) 


10 69 +lul 


8 - 64 -L56 


0.40^* 


? 14+1.40 


47.05 


< 1.01 



a The rest-frame Lya EWs are determined by assuming 
Values in parentheses indicate the EW in case of using 
Lya wavelength to estimate the continuum. 
b The dumpiness parameter (see text for further details). 



that all LAEs are at z = 4.86 except for #35915. 
broadband images both redward and blueward of 
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Figure 9. Observed SEDs and the best-fitting models for group I — III LAEs. 
In each panel, the observed SED is shown with filled circles; the best-fitting 
model SED is indicated by opened squares. Solid line represents the best- 
fitting model spectrum with Ha emission line. Arrows indicate the 2tr upper 
limits. V-band photometry is also shown in the figures but is not used in SED 
fitting. 



ject. For LAEs with upper limit, we used the upper limit to 
estimate the magnitude. Despite the small size of the sample, 
the rest-frame optical luminosity is likely to be a tracer of the 
stellar mass. The relationship between the stellar mass and 
extinction-corrected optical absolute magnitude also shows a 
correlation (Figure [TOl b)). A mass-to-light ratio is not per- 
fectly linear; it becomes larger in the brighter LAEs. This 
trend is similar to the results of LBGs at the same redshift 
bv lYabe et al. I (|2009) which will be discussed in more details 
in the next section. From the relationship between the stellar 
mass and the rest-frame optical luminosity, it is implied that 
most LAEs at z = 4.86, which could not be detected above 
2<T limiting magnitudes in 4.5^m band (TableHJ, are likely to 
have stellar masses lower than 10 9 M Q . 

The best-fitting stellar ages are in a wide range from 7.4 
Myr to 437 Myr with the median age of 25 Myr. Although 
we did not limit the model age at the age of the universe at 
z = 4.86 (~ 1 .2 Gyr), there is no LAE with the best-fitting age 
exceeding the cosmic age. The typical error on age of each 
object is very large (~ 1.7 dex). The uncertainties of ages 
cover most of the age range available in the models. Thus it 
is difficult to draw any conclusion on the correlation of them. 

The color excesses E(B-V) range from 0.14 mag to the 
value as high as 0.40 mag. The median value is 0.27 mag. 
The typical error on E(B-V) is ~ 0.2 mag. It is seen in 



Figures \TU\ c) and (d) that dust extinction seems to show 
some correlation with both rest-frame UV and optical mag- 
nitudes. However, the trends are less significant if we con- 
sider the uncertainties and upper limits. Similarly, no trend 
between the color excess and other fitting properties can be 
seen when we take the uncertainties into account (Figure ITTb. 
As mentioned in the last part of section [3] IRAC photometry 
of group III LAE (#23098) may be contaminated. By remov- 
ing the LAE from group III, the derived color excesses are 
0.14-0.30 mag; any correlations between the dust extinction 
and magnitudes or other properties cannot be seen. The pres- 
ence of some amount of color excess indicates that LAEs at 
this redshift are not free from dust, whic h is consistent with 
what lFinkelstein et al. I d2007ll2008l 1200 9) found for LAEs at 
z ~ 4.5. Many a uthors (Neufeld 1991; HanieiL&Oh] 2006; 
Finkelstein et al. 2009) proposed that the presence of dust in 
clumpy clouds can enhance the Lya EW of a galaxy. To test 
this hypothesis, we estimated the dumpiness parameter (q): 
qr is Lya line opac ity where r is an optical depth for the 
continuum (Finkelste in et al. 1 120081) . If q < 1, the Lya EW 
is enhanced on the hypothesis that Lya photons suffer less 
extinction than the continuum photons and the clumpy cloud 
model would be supported. We calculate q from the ratio 
of observed and intrinsic Lya luminosities. The observed 
Lya luminosity is derived from the observed Lya EW and 
the continuum flux density, while the intrinsic one is calcu- 
lated from the Ha luminosity by assuming case B recombi- 
nation. Except for #12766 and #23098 that show upper limits 
on the q parameter, all LAEs show q > 1 (Table |2]l suggest- 
ing that Lya photons suffer from dust extinction larger than 
the continuum. The uncertainty of q value is large; neverthe- 
less, the q values are all near 1-2 suggesting that although 
the Lya emission appears to be more attenuated than the con- 
tinuum, the difference is not so much. In section [3] however, 
we mentioned that the observed Lya EWs may be underesti- 
mated if the Lya emission line does not fall into the central 
part of NB1\ 1 filter and are further underestimated due to the 
non-detection in I c band for some LAEs. Underestimation of 
the Lya EW results in the overestimation of the q parameter. 
Accordingly, the objects without spectroscopic confirmation, 
i.e., #30702 and #22962, could have q values less than those 
shown in Table|2] if their Lya line did not fall into the central 
part of the NB1 1 1 filter. The spectroscopic observations are 
desirable to investigate the real Lya EWs and q values of the 
LAEs. 
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Figure 10. Plots of derived stellar masses and color excesses of group I — III LAEs at z = 4.86 against the rest-frame UV or optical absolute magnitudes. Filled 
circles, open triangles, and open circles represent group I, group II, and group III LAEs, respectively. Vertical error bars represent the 68% errors as described in 
text. Arrows indicate the objects whose absolute magnitudes were calculated from the upper limits. 



The best-fitting star formation rates (SFRs) of LAEs range 
from 55 M Q yr _1 to 209 M Q yr _I . The median value of the 
derived SFRs is 132 M Q yr _I . The typical error on SFR for 
each object is ~ 0.9 dex. These SFRs are much higher than 
those of LAEs from other studies (e.g., Gawiser et al. Il2007t 
iNilsson et al. 1 12007). This is probably due to the selection ef- 
fect. Our SED fitting sample has to be bright enough to be 
detected in either ground-based broadbands or IRAC bands 
or both. Thus the selected LAEs are more luminous which 
results in the higher SFRs. The relation between the SFR 
and other fitting properties are shown in Figure [TTT c). (e), and 
(f). Looking only at the best-fitting values, we see that SFR 
seems to show correlation with stellar mass and dust extinc- 
tion. However, the correlations are not significant when we 
consider the uncertainties or exclude group III LAE. 

Since some of the LAEs show very young ages, an effect of 
nebular continuum emission may not be negligible. We exam- 
ine the effects of nebular emission on the SED fitting results 
and describe more details in Appendix B. Briefly, we found 
that stellar mass estimation of all LAEs is robust (±0.04 dex). 
Age, color excess, and SFR for group I LAEs are also ro- 
bust (±0.24 dex, ±0.02 mag, and ±0.10 dex, respectively). 
For group II and III LAEs, age, color, and SFR do not agree 
well with each other, but agree within their large uncertainties. 
Ha emission accounts for the largest expected nebular line 
emission to be detected from our LAEs. Including Ha emis- 
sion to the model spectrum when we performed SED fitting in 



the last section is likely to be the reason why adding the other 
nebular emissions give the similar results. 

7. COMPARISONS TO LBGS AT Z ~ 5 

In this paper, we aim at comparing the derived stellar pop- 
ulations of LAEs at z ~ 5 to those of LBGs at t h e same red- 
shift and in the same field bv lYabeetaLl (12009ft . lYabeetal.l 
(2009) used LBGs at z ~ 5 with z' < 26.5 mag selected by 
llwata et al~ (2007); they used the same data set and the same 
V-dropout criteria as those used in this study. In addition, the 
SED models and the fitting technique are also the same. In 
order to make a fair comparison, we used only LAEs whose 
z 'totai magnitudes are brighter than 26.5 mag. Because part of 
the LBGs are spectroscopically confirmed with the mean red- 
shift of z~ 4.7 (Ando et al. 2007; Kajino et al. 2009), taking 
the same limit on z' magnitude is equivalent to taking the same 
limit on rest-frame UV luminosity. There are six LAEs with 
z 'totai < 26.5 mag in the IRAC isolated sample. Two of them 
are the group I LAEs; one is in group II (#12766). The oth- 
ers are group IV LAEs for which we could not perform SED 
fitting. Therefore, we used three LAEs from group I and II 
for the comparis on. Among t hese three LAEs, one (# 35915) 
appear s in llwataet al. I ([2007) catalog and also in Yabe et al. 
(2009) . Though w e used the same set of broadband data as 
llwata et al71 d2007l) to select the sample, the images in our 
case were smoothed to have the same seeing size as that o f 
the NB1W image before making photometry (section 12.1.21 ). 
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Figure 11. Relations of output parameters from SED fitting for group I — III LAEs at z = 4.86. Filled circles, open triangles, and open circles represent group I, II 
and III LAEs, respectively. 



Accordingly, magnitudes measured from the images before 

and after smoothing are not n ecessarily the same . 

From 617 LBGs at z ~ 5, dlwata et al. 1120071) . lYabe et al. I 
(2009) selected 170 LBGs by eye inspection as isolated ob- 
jects both in z'-band and in IRAC images. The 170 LBGs 
are divided into four categories according to the detection in 
3.6/im and 4.5/im filters. The LBGs that are detected in both 
3.6/im and 4.5/im images are in category 1, while those only 
detected in either 3.6/im or 4.5/im are put in category 2 or 3, 
respectively. The rest that are detected neither in 3.6/im nor 
4.5/im images are in catego ry 4. We used categories 1, 2, and 
3 LBGs in the comparison. lYabe et aLl (2009) performed the 
SED fitting for 64 LBGs in category 1 . For categories 2 and 
3 LBGs, we refitted them by using an upper limit either in 
4.5/im or 3.6/im band, respectively. 17 LBGs are best fitted 



with the stellar ages larger than the cosmic age at z ~ 5 (~ 1.2 
Gyr). The reason for the large ages is considered to be the 
assumption of constant star formation history or low-z inter- 
lopers (Yabe et al. 2009; Sawicki & Yee 1998; Papovich et 
al. 2001). The derived ages are generally older than those 
derived by assuming other star formation histories. We thus 
excluded these objects in the following comparison. We even- 
tually have 88 LBGs from categories 1, 2, and 3 to compare 
with our three LAEs from groups I and II. It should be noted 
that since the NB7 1 1 filter covers only a small fraction of the 
redshift range, LAEs that are not in this redshift range are ex- 
pected among the LBG sample. 

Before going to the comparison between stellar populations 
of LAEs and those of LBGs, we try to figure out what their 
photometry can tell us without SED fitting. Another thing 
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indicated by red filled circles and red open triangle, respectively. LBGs in category 1 are plotted in filled triangles; those in categories 2 and 3 are plotted with 
black open triangles. Arrows show the limits in magnitudes or colors. Errors of LBGs are taken from Yabe et al. 12009). 



which should be kept in mind is that three LAEs which are 
used to compare with the LBGs in this section (i.e., group 
I LAEs and #12766, hereafter "compared LAE sample") are 
not necessarily good representatives for all LAEs at z = 4.86. 
As shown in Figure [8] the compared LAE sample (grey his- 
togram) is relatively bright and has comparatively red UV 
color as compared to the whole LAE sample. Figure Q~2] 
shows the color-magnitude diagrams of the compared LAE 
sample and the LBGs at z ~ 5. No clear trend is seen be- 
tween z' magnitude and I c -z' color (Figure[T2Ta)). The figure 
shows that the LAEs occupy the faint part of z'-magnitude 
distribution and the red part of I c -z' color distribution of the 
LBGs. Distribution of these LAEs at the faint part of UV lu- 
minosity among LBGs may support the deficien cy of stron 
Lya emission among bright LBGs claimed by Ando 



strong 
et al. I 



(2006), though the sample size is only three. Figure [T2T b) 
shows the color-magnitude diagram of z' - [4.5/im] color and 
4.5/im magnitude. The z 1 - [4.5/im] colors of LBGs spread 
over a wide range from -1 to 5 mag. The LAEs are on the 
distributions of LBGs, but they locate in the slightly fainter 
part in 4.5/im-band magnitude and the relatively bluer part of 
the color distributions than the LBGs. 

Output parameters of SED fitting are plotted against the 
rest-frame UV and optical luminosities in Figure Q~3] The 
top panels show the comparisons between the derived stel- 
lar masses of LAEs and those of LBGs. The stellar masses 
of the compared LAE sample range from 10 9 M© to 10 I0 M Q . 
In contrast, LBGs distribute in a larger range from 1O 8 M to 
1O H M0. It is seen in the figure that the LAEs locate in the 
region where LBGs distribute and in the almost middle part 
of LBGs' distribution. At the same bin of rest-frame UV or 
optical luminosity, LAEs seem to have comparable masses to 
LBGs. 

Showing no relation between ages and luminosities, Figures 
Qjjc) andQjfd) indicate that the LAEs again lie on the same 
distribution of LBGs. The median ages are similar; they are 
25 Myr and 22 Myr for LAEs and LBGs, respectively. How- 
ever, age uncertainties of the LAEs cover the whole range of 
age distribution of the LBGs. It is difficult to draw any con- 
clusion. 

Figures [I"3f e) and fOl f) show the plots of the derived color 
excesses against UV and optical absolute magnitudes, respec- 
tively. Although the dust extinction of the LAEs is not zero, 



they seem to lie at the relatively lower region of the color ex- 
cess distribution of LBGs. At the fixed rest-frame UV or op- 
tical luminosity, there seems to be no difference in E(B - V), 
even if we take the uncertainties into account. Figure [Di e) 
shows that there is no relation between the dust extinction and 
the rest-frame UV luminosity. In contrast, Figure [T3l f) seems 
to show a correlation of LBGs between the color excesses and 
the optical magnitudes. Though we can not state any correla- 
tion for only 3 LAEs, the LAEs still lie on the correlation of 
LBGs. 

The SFRs are plotted against rest-frame UV and optical ab- 
solute magnitudes in Figures [T3l g) and [L37 h). respectively. 
Both figures show a correlation of SFRs of the LBGs with 
UV (but weak) and optical luminosity. It is seen from the fig- 
ures that the LAEs locate in the lower part of the LBGs' dis- 
tribution. A median SFR of the LAEs is 132 M Q yr _1 , while 
it is 187 M yr _1 for LBGs. Note that the LAEs used in the 
comparison are biased toward the brightest sample among the 
whole LAE sample (section|4|i, which probably results in the 
higher SFRs of the LAEs. In conclusion, no significant dif- 
ferences in the stellar properties between LAEs and LBGs are 
seen at the same luminosity. 

A physical reason that makes an LAE an emitter is still 
unclear. In order to investigate whether or not the dif- 
ference between LBGs and LAEs has any dependence on 
the equivalent widths (EWs) of Lya emission, we plotted 
the output parameters of SED fitting against the rest-frame 
Lya EWs in Figure[l4] Group I— III LAEs a re plotted with five 
LBGs which have spectroscopic Lya EWs dAndo et al. I2004t 
Kaiino et al. 2009). Arrows in the figures represent the lower 
limits on EWs of the LAEs; I c magnitudes of them are fainter 
than 2cr limiting magnitude. Note that, except for one spectro- 
scopically confirmed LAE, the Lya EWs of LAEs are likely 
to be lower limits if their Lya emission lines do not fall into 
the center of NB711 band. According to the figure, we can- 
not find any significant correlation between the stellar prop- 
erties and the rest-frame Lya EWs. Recently, Kornei et al. 
(2009) studied these relations for z ~ 3 LBGs with spectra. 
They found that the stellar mass does not correlate with the 
Lya EW, while large Lya EW is seen in older, lower SFR, 
and less dusty LBGs. Such trends are difficult to isolate with 
only three LAEs. A larger sample at z ~ 5 is required to see 
clearer relations if any and to be definitive. 
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8. SUMMARY 

In this paper, we studied the stellar properties of Lyman 
alpha emitters (LAEs) at z = 4.86 using SED fitting. By 
narrowband and broadband observations by Suprime-Cam 
on Subaru telescope, 24 LAEs were selected in the area of 
~ 508.5 arcmin 2 around GOODS-N field. In addition to the 
optical photometry, we obtained the mid-infrared photome- 
try from data taken by IRAC on Spitzer space telescope in 
the GOODS-N field as well as the surrounding area in or- 
der to cover most part of Subaru area. We selected 12 LAEs 
that are isolated from the neighboring objects. We performed 
SED fitting of five LAEs which are detected above 2a mag- 
nitude limits in more than 2 bands. Selecting those five LAEs 
could introduce a bias toward the bright red galaxies. Model 
SEDs are built by assuming the constant star formation his- 
tory with fixed metallicity of 0.2Z Q , the Salpeter IMF ranging 
from 0. IMq to 100 M Q , and the extinction law of Calzetti et 
al. (2000). The derived stellar masses of the LAEs range from 
10 8 to 10 10 M Q with the median value of 2.5 x 1O 9 M . The 
derived ages cover wide range from 7.4 Myr to 437 Myr with 
the median value of 25 Myr. The color excess are between 
0.1 -0.4 mag, indicative of the presence of some amount of 
dust. Star formation rates (SFRs) are in the range of 55 - 209 
M Q yr _1 . The median color excess and SFRs are 0.27 mag and 
132M Q yr _1 , respectively. The high SFRs are probably due 
to the selection effect; we selected the LAEs that are bright 
enough to be detected in rest-frame UV and optical bands, 



which results in selecting the LAEs with the higher SFRs. We 
investigate the correlations between the stellar properties de- 
rived by SED fitting and the photometric properties of LAEs 
and found no significant correlation due to both small size of 
the sample and the large uncertainty on fitting results. 

The main objective of this study is to compare LAEs to 
LBGs at the same redshift. LBGs were selected by V-dropout 
criteria (Iwata et al.2007). Their stellar populations were de- 
rived by Yabe et al. (2009). Because those LBGs are se- 
lected from the same set of data and stellar population was 
derived by the same SED fitting model, we can make a fair 
comparison between their stellar populations. We compared 
three LAEs to 88 LBGs down to the same UV luminosity 
limit. These three LAEs are the brightest and reddest ones 
among the whole LAE sample in this study. The comparisons 
of SED-fitting parameters show that LAEs locate in the re- 
gion where LBGs distribute; the physical properties of LAEs 
and LBGs occupy similar parameter spaces. At the same rest- 
frame UV or optical luminosity, there is no difference of stel- 
lar properties between LAEs and LBGs. In order to figure 
out properties that control the Lya EWs, we plotted the out- 
put parameters against the rest-frame Lya EWs ranging from 
to 40 A. We could not find any significant correlation be- 
tween them. A larger number of sample is needed to see any 
correlation of Lya EW if exists. 
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APPENDIX 
EFFECTS OF NEBULAR EMISSION 



In order to investigate the effects of nebular emission on the synthesis models, we used the stellar population synthesis model by 
Fioc & Rocca-Volmerange (1997, 1999; hereafter PEGASE). The PEGASE can add nebular emission (continuum and emission 
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Table 3 

The best fit results for PEGASE model 



ID 


Field 


log[Mass] 
(Mq) 


log [Age] 

(yr) 


E(B-V) 

(mag) 


log[SFR] 
(Mgyr 1 ) 


x 2 . 


Group I 


30702 
35915 


GOODS-N 
GOODS-N 


q 49+1.13 

-Q-Z2 

y - zy -0 17 


7 J 2+2.84 

7 90+2:06 


30 +,u2 
-Q-3D 

20+ 8 ' 2 


2 32 + ^'^ 

2 03+^ 
z - UJ -i 10 


0.001 
0.05 


Group II 


12766 
22962 


GOODS-FF 
GOODS-N 


9 25+ 1M 
853+^ 


6 34+ J% 
-0-34 

6 95+ 2 - 4 ' 


Q3J+U.17 

-Q-3I 

13+ 8 " 


2 90+" ' ' 

1 59+ T: ® 
'■"-1 m 


44.88 
46.64 


Group III 


23098 


GOODS-FF 


10.29^ 


1 38+2 L .>2 


56 +u - J 
u - JD -0.44 


9 Q4.+ 1^3 
z - yH -1.78 


47.05 



lines) to the model spectrum. The hydrogen emission lines are computed from the number of Lyman continuum photons by 
assuming case B reco mbination. Other emission lines are calculated from the observed ratios to H/3 for typical local starbursts 
(Fioc & Rocca-Volmerange 1997). Other prescriptions are the same as those used in the case of BC03 model described in the 
main text. We used Salpeter IMF (1995) with lower and upper mass cutoffs of 0.1 and 100 Mq, respectively, adopted metallicity 
of 0.2Z Q , and assumed the constant star formation history. We re-scaled the age step of PEGASE into the logarithmic scale as we 
did for BC03 model. The effect of dust attenuation is taken into account by using Calzetti extinction law (Calzetti et al. 2000). 
E(B-V) varies from 0.0 mag to 0.8 mag by 0.01 mag step. The attenuation due to IGM is followed by Madau (1995). The 
redshift is fixed at z = 4.86 except for #35915 which has a spectroscopic redshift (section[5]l. 
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Figure 15. Observed SEDs and the best-fitting models for LAEs in group I (left panel), group II (middle panel), and group III (right panel). In each panel, the 
observed SED is shown with filled circles; the best-fitting model SEDs are indicated by opened squares and asterisks for BC03 and PEGASE models, respectively. 
Black line represents the best-fitting BC03 model spectrum with Ha emission line, while blue one represents the PEGASE model. 2cr upper limit of V-band 
photometry is also shown with the arrow; however, the data is not used in SED fitting. 



The results obtained for group I— III LAEs with PEGASE are summarized in table [3] The best-fit spectra are shown as a 
comparison with BC03 best-fitting spectra in Figure [15] The comparisons of resulting parameters are presented in Figure [16] 
Figure [ToTa) shows that the stellar masses derived by both models are well in agreement with each other. The stellar masses 
obtained by the PEGASE model are averagely smaller than those by the BC03 model by about 0.04 dex. As seen in Figure fTST b). 
the stellar ages of most LAEs derived by PEGASE are younger than those derived from BC03 models. Inclusion of nebular 
emission averagely decreases the ages of all LAEs by 0.46 dex. However, due to the large uncertainties in determining the stellar 
ages, the derived ages seem to be in agreement within the uncertainty range. Figure [TSTc) shows that the color excesses derived 
by PEGASE are comparable to those by BC03 model. The average difference is 0.04 mag except for the group III LAE. As 
seen in Figure [TBT d). the SFRs derived by both models seem to be comparable for most of the LAEs. However, the rather large 
differences are seen for the group III LAE and one in group II (#12766). SFRs increase averagely about 0.28 dex when nebular 
emission is included. 



16 



Yuma et al. 



13 



10 

CD 
W 

S10 12 

as 
o. 

^10 11 

CO 



10 



10 



CO 

E 
_cfl 
"53 
CO 



10 s 
10 8 
10 7 



r (a) ■ 


- 




/ 


r 




/ - 








A 


Y , 


■ 














- 


/ 




_ 


■/ 



10' 10 c 



10 9 10 10 10 11 10 12 10 13 




Stellar mass (Msun) BC03 



10° 10 
Stellar age (yr) BC03 



10 10 10 11 



CO 
O) 
CD 

o_ 

> 

CO 



1 

0.8 
0.6 
0.4 
0.2 




- (c) 


i ' i ' / 


- 




i , i 










C — i 



0.2 0.4 0.6 0.8 
E(B-V) BC03 



10 H 



CD 
W 

to 1f1 3 
O) 10 

CD 

CL 



10" 



C10 1 

GO 



10 u 



3 


i /'- 




- 


"'/' • 


- i 


• 


-I -, 
I 


r / 


1 1 

Mill 



10 u 



10' 



10" 



10 J 



SFR (Msun/yr) BC03 



10" 



Figure 16. Comparison of output parameters from BC03 and Pegase models. Filled circles, open triangles, and open circles represent group I, group II, group 
III LAEs, respectively. 



